Electron cooling and Debye- Waller effect in photoexcited bismuth 
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By means of first principles calculations, we computed the effective electron-phonon coupling 
constant Go governing the electron cooling in photoexcited bismuth. Go strongly increases as a 
function of electron temperature, which can be traced back to the semi-metallic nature of bismuth. 
We also used a thermodynamical model to compute the time evolution of both electron and lattice 
temperatures following laser excitation. Thereby, we simulated the time evolution of (1 -1 0), (-2 
1 1) and (2 -2 0) Bragg peak intensities measured by Sciaini et al [Nature 458, 56 (2009)] in 
femtosecond electron diffraction experiments. The effect of the electron temperature on the Debye- 
Waller factors through the softening of all optical modes across the whole Brillouin zone turns out 
to be crucial to reproduce the time evolution of these Bragg peak intensities. 

PACS numbers: 63.20.K-, 61.80.Ba, 78.47J- 



Femtosecond time resolved X-ray or electron diffrac- 
tion emerges as a unique tool to study laser-induced 
structural dynamics with a picometer spatial and fem- 
tosecond temporal resolution. A large class of systems 
including organic compounds [![, correlated materials 0, 
3] and semi-metallic systems like graphite]!, Q or 
bismuth[6l-lll| have been investigated, providing enlight- 
ening information about atomic motions following laser 
excitation. However, it is often difficult to interpret the 
time evolution of all Bragg peak intensities and to un- 
derstand the different time scales showing up in these 
experiments because models or calculations describing 
laser-induced electron-phonon processes are lacking. 

A first step towards a better understanding of time 
resolved X-ray diffraction experiments performed on 
bismuthQ was recently achieved[l2j]. Indeed, a thermo- 
dynamical model based on the assumption that an elec- 
tron temperature T e and a lattice temperature T; can be 
defined at any time delay following the laser pulse arrival, 
was shown to reproduce the (111) bragg peak intensities 
measured by Fritz et aZ0] for time delays up to ~ 2.5 ps. 
While most of the parameters encompassed in the model 
were obtained from first principles calculations, the ef- 
fective electron-phonon coupling constant Go governing 
the electron cooling was taken as a temperature indepen- 
dent parameter fitted to reproduce the measured Bragg 
intensities. Obviously, knowledge of the temperature de- 
pendence of Go is mandatory to perform simulations over 
longer time scales. 

In this letter, we computed the phonon spectrum of Bi 
for a large range of electron temperatures and calculated 
the effective electron-phonon coupling constant Go(T e ) 
using the theory developped by Allen[13] which is applied 
for the first time to a material showing a strongly varying 
electron density of states at the Fermi level. We found 
that all the optical modes are affected by an increase of 



T e and that Go strongly increases as T e increases. Then, 
we solved the model proposed in Ref. [12| and simulated 
the intensity decay of Bragg peaks measured by Sciaini 
et a/pjj in femtosecond electron diffraction (FED) ex- 
periments for time delays up to 14 ps. Our calculations 
firmly establish that the softening of all optical modes 
is crucial to achieve a good agreement between experi- 
ment and theory. This finding parallels that of Johnson 
et «?[Io| who scale all phonon modes by an empirical 
constant factor smaller than one to reproduce their nor- 
malized diffracted intensities. 



All the ab-initio calculations were performed within 
the framework of the local density approximation (LDA) 
using the ABINIT code[I3. s P in 

-orbit coupling was 
included and an energy cut-off of 15 Hartree in the 
planewave expansion of wavefunctions as well as a 16 x 
16 x 16 kpoint grid for the Brillouin zone integration 
were used. We computed the evolution of the phonon 
spectrum as a function of T e using density functional 
perturbation theory (DFPT)[15]. Our calculations show 
that the phonon spectrum is not affected by T e below 
600 K. However, an increase of T e above this critical 
temperature leads to a redshift of all optical phonon 
modes. Figure 1(a) shows the calculated phonon density 
of states for T e = K and T e = 2100 K. All the optical 
modes are softened for T e = 2100 K whereas the acoustic 
modes are practically unaffected. This effect has already 
been highlighted in other theoretical works [H, Oj} and 
should be contrasted with the bond hardening predicted 
theoretically [l8[ and observed experimentally in gold[l9j]. 



Knowledge of the phonon spectrum is essential to com- 
pute the decay of electron energy due to electron-phonon 
interaction. A key dimensionless quantity is the so called 
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generalized Eliashberg function 13] defined by 
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where the sum is performed on ATj. (iV 9 ) electron 
(phonon) wavevectors over the Brillouin zone, g{t F ) is 
the density of states per unit cell at the Fermi level, ekn 
are Kohn-Sham energies for band n and wavevector k 
and where the electron-phonon matrix elements read 
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where (q, A) is the displacement of the p th atom of 
mass M in the direction a for the mode (q, A) of fre- 
quency ujq\. By solving the Boltzmann equation for both 
the electron and phonon systems, one can show that the 
variation of electron energy per unit volume is given by 



dt 



G (T e ) x (T, - T e 



(3) 



provided that the lattice temperature TJ is larger than 
the Debye temperature Op ~ 119 K of Bi. The effective 
electron-phonon coupling constant Go(T e ) reads [l3T| 



G (T e ) 
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where i> is the unit cell volume, fc^ is the Boltzmann con- 
stant, and /(e) is the Fermi-Dirac distribution at tem- 
perature T e . This expression shows that the integra- 
tion over e should be done in a range of energy around 
e F whose width increases with T e . The usual approx- 
imation is tantamount to neglecting the energy depen- 
dence of a 2 F(e,e + huj,uj) which is then approximated 
by a 2 F(e F , e_F, w) = a 2 F(ui) when the phonon energies 
are also neglected. Thus, the electron-phonon coupling 
constant becomes independent of T e and is given by 



G = —k B X(uj 2 )g(e F ), 
vn 



(5) 



where A(w 2 } is the second moment of the Eliashberg func- 
tion a 2 F{uj). Eq. [5]is justified for low electron tempera- 
tures or for all electron temperatures when the electron 
density of states (DOS) is only weakly energy dependent. 
The last assumption is far from being satisfied since Bi 
is a semi-metallic material with a strongly varying DOS 
near the Fermi level. Therefore, a temperature depen- 
dent Go is expected. 

The ab-initio calculated electron-phonon coupling con- 
stant Go is shown in Fig. [U(b). The low temperature 
value is well described by Eq. [S] which explains why Go 
is so small since g(e F ) ~ 1.94 x 10" 2 eV _1 per cell. The 
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FIG. 1: (color online) a/ Calculated phonon density of states 
in meV" 1 per cell for T e = K (thin line) and T e = 2100 
K (thick line). The inset shows the fraction of acoustic 
(TA+LA) and optical (TO and LO) phonons contributing to 
the effective electron-phonon coupling Go for electron temper- 
ature T e ranging from 294 K to 4000 K. b/ Calculated effective 
electron-phonon coupling constant Go (in W.m _3 .K _1 ) as a 
function of T e (in K) compared to the values fitted in Ref.[l2T] 
to reproduce the time evolution of (111) Bragg peak intensi- 
ties measured by Fritz et ai[3 • The inset shows the calculated 
bandstructure of Bi along high symmetry lines for T e = K. 



huge increase of Go as T e increases can be easily inter- 
preted by looking at the electronic structure of Bi de- 
picted in the inset of Fig. [TJb). At low T e , only intraband 
scattering processes involving nearly zone center phonons 
or interband scattering processes involving phonons con- 
necting the L pocket to the T pockets participate to the 
energy exchange between the electrons and the lattice. 
Thus, the energy is transferred to only a few phonon 
modes in this regime. When T e increases, the number 
of phonon modes participating in the intraband scatter- 
ing processes increases tremendously because of larger 
q-wavectors allowed by momentum and energy conserva- 
tion. Consequently, Go increases by four orders of magni- 
tude when T e increases from 100 K to 4000 K. The inset 
in Fig. QJa) shows that both optical and acoustic phonons 
contribute to electron cooling. At 2100 K, transverse op- 
tical phonons (TO), longitudinal optical phonons (LO) 
and acoustic phonons (TA+LA) respectively contribute 
for 51.6, 27.1 and 21.3 % to the total electron-phonon 
coupling constant Go- In Ref.[l2j], the electron-phonon 
coupling constants were fitted to reproduce the time evo- 
lution of (111) Bragg peak intensities measured by Fritz 
et a/[7j for four different pump laser fiuences. All the 
fitted values denoted as stars in Fig. [TJb) are underesti- 
mated by ~ 20% with respect to the calculated values. 
However, the (111) Bragg peak intensities are still per- 
fectly reproduced by using our calculated Go(T e ), which 
was also succesfully used to reproduce the time evolution 
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FIG. 2: (color online) a/Mean square displacements (in A 2 ) 
of one Bi atom as a function of lattice temperature T; (in 
K) from ab initio calculations for T e =294 K and T e =2100 K 
(thin and thick solid lines). The vertical line indicates the 
melting temperature of Bi. b/Calculated normalized intensi- 
ties defined by I hkl (T h T e = T ; )/I hk i(T ; = 294 K,T e = 294 K) 
for (1 -1 0), (-2 1 1) and (2 -2 0) Bragg peaks as a function 
of lattice temperature Ti in K. The horizontal lines are the 
experimental intensities obtained in FED experimentsfrO] 14 
ps after the arrival of the 200-fs pump pulse of fiuence 0.84 
mj.cm -2 . The experimental data are shown in Fig. EJc). 



of the electron temperature inferred from time-resolved 
photoemission experiments [2fj| . 

We next calculate the mean square displacement of 
atoms with respect to their equilibrium positions as a 
function of both lattice temperature 7] and electron tem- 
perature T e . Defining the displacement of the p*' 1 atom 
in the unit cell in the direction a as ajj, one can show 
that 



q,A 1 
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where riq.\ is the Bose occupation factor at temperature 
Ti for a phonon with frequency w q A- The mean values 
defined by Eq. [5] have been calculated using the method 
introduced by Lee and Gonze[2l|- Our calculations show 
that the anisotropy in the mean square displacements 
is negligible. Figure |2Ja) displays the mean square dis- 
placement of one Bismuth atom (in A 2 ) for T e = 294 
K as a function of Tj (thin solid line). The root mean 
square displacement increases from 0.070 A at T z =0K 
to 0.241 A at T/=294 K. The room temperature value is 
slightly overestimated with respect to the experimental 
value of 0.21 A extracted from LEED experiments 22 1 . 
The softening of the optical phonon modes over the whole 
Brillouin zone sketched in FigfTJa) has a strong impact 
on the mean square displacements which are as higher 
as the phonon frequencies are lower (see Eq. [5]). Figure 



[2ja) shows the increase of the mean square displacements 
upon increasing T e from 294 K (thin solid line) to 2100 
K (thick solid line) and suggests that the Bragg peak 
intensities are affected by T e . 

The diffracted intensity for a given scattering vector 
defined by the Miller indices (h,k,l) is given by 



Ih,k,i(Ti, T e ) = J? fc ,exp[-2W(T,,T e )] 



(7) 



where the intensity I^ki' which can be affected by the 
coherent A lg phonon coordinate [HI, is reduced by the 
so-called Debye- Waller factor. Here W(Ti,T e ) reads 



W(T h T e ) 



a.B 



(8) 



where 



is obtained from equation [6] Figure [2fb) 
displays the calculated ratios Ihkl(7j,T e = T/)/Ihki(T] = 
294 K, T e = 294 K) for (1 -1 0), (-2 1 1) and (2 -2 0) Bragg 
peaks as a function of T; . The horizontal lines correspond 
to the normalized intensities obtained by Sciaini et a/[ll[ 
in their FED studies 14 ps after the arrival of the laser 
pulse. As can be seen from Fig|3]Jc), the experimental 
normalized diffracted peak intensities are stationary for 
this time delay. Assuming that the lattice is at equi- 
librium with the electron system for larger time delays, 
we can deduce from our calculated normalized intensities 
three different lattice temperatures shown by arrows in 
FigJ^b) which are close to each other. The temperature 
Ti >eq reached by the lattice, when averaging these three 
values, is found to be 408 K and is smaller than the tem- 
perature T^ eq ~ 460 K obtained by Sciaini et alYLM on 
the basis of a parametrized Debye- Waller model f23j. 

In order to simulate the time-resolved FED experi- 
ments performed by Sciaini et ctZflll] on a free standing 
Bi film of thickness L=30 nm excited by a near-infrared 
laser pulse whose FWHM is t w = 200 fs, we used the 
model proposed in Rcf.^ljl where the source term P(z, t) 
describing the energy deposited per unit volume and time 
by the laser pulse at time t and at position z reads 



P(z,t) 
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where Fi nc is the incident fiuence and where g(z) is a 
function plotted in the inset of Fig(3Ja) and calculated 
by using both the imaginary and real part of the op- 
tical index of Bi for a wavelength A = 775 nm. We 
note in passing that the absorption coefficient A de- 
fined by Jq g(z)dz is slightly larger for a 30 nm Bi film 
(A = 0.36) than for a semi-infinite film (A = 0.31). 
The incident fiuence Fi nc needed to reproduce the lat- 
tice temperature rise extracted from Fig. [2Jb) is given 
by F mc = L x Q x {T leq - T lfi )/A where T lfi is the ini- 
tial lattice temperature and C\ is the Dulong-Petit lattice 
specific heat. Thus, we estimate an incident fiuence of 
1.11 mJ.cm~ 2 . This value is overestimated by 27 % with 
respect to the measured incident fiuence of 0.84 mj.cm , 
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FIG. 3: (color online) a/ and b/ Lattice Tj and electron tem- 
perature T e in K, as a function of time delay t in ps obtained 
from the model proposed in Ref . [12| for a 30-nm-thick Bi film 
excited by a 200-fs laser pulse of incident fluence Fi„ c =l.ll 
mj.cm -2 at a wavelength of 775 nm. The inset shows the 
function g in nm -1 as a function of z in nm. c/ Normal- 
ized intensity decays of (1 -1 0) , (-2 1 1) and (2 -2 0) Bragg 
peaks measured by Sciaini et aZ[l]| as a function of t com- 
pared to the theoretical calculations. The thin lines repre- 
sent I h , k ,i(T h T e = Ti)/I h , k ,i(Ti = 294 K,T e = 294 K) while 
the thick lines represent I h , k ,i{ T i-. T e) I 'h,k,l(Ti = 294 K > r * = 
294 K). 



which might be traced back to the problem of making an 
accurate measurement of the laser fluence. 

Using our calculated effective electron-phonon cou- 
pling constant Go shown in Fig. QJb), we numerically 
solved the three coupled differential equations underlying 
the model proposed in Ref. [12j . Thereby, we obtained the 
spatial and temporal evolution of the coherent phonon 
coordinate u and of both electron T e and lattice temper- 
ature Ti following the arrival of the 200 fs laser pulse on 
the sample. As shown in Fig. EJb), the electron temper- 



overheating of the electron system, the electron tempera- 
ture decreases whereas the lattice temperature increases. 
At t = 14 ps, the lattice is not at equilibrium with the 
electron system due to the slow down of electron cool- 
ing for low electron temperatures. However, the lattice 



temperature almost reaches T) 



408 K because the 



ature reaches its maximum value T e 



2257 K only 



0.3 ps after the arrival of the laser pulse. After a strong 



energy stored in the electron system is very small. 

The normalized intensities of (1 -1 0), (-2 1 1) and (2 
-2 0) bragg peaks measured by Sciaini et al [ll[ in FED 
experiments are shown in Fig. EJc). These intensities are 
compared to the theoretical intensities (thick solid lines) 
defined by I h , k ,l( T h T e)/Ih,k,l( T l = 294 K, T e = 294 K). 
The agreement between theory and experiment is note- 
worthy given the fact that all parameters of the model 
have been obtained from ab-initio calculations. In order 
to point out the role played by T e , we have also displayed 
Ih,k,i{TuT e = Ti)/I h>k ^Ti = 294 K, T e = 294 K) (thin 
solid lines) in Fig. EJc). Our results show that it is not 
possible to reproduce the normalized intensity decays if 
we assume that the electron system is at equilibrium with 
the lattice. The Bragg intensities are systematically over- 
estimated with respect to experiments. From a physical 
point of view, the softening of the optical phonon modes 
across the whole brillouin zone when T e is larger than 600 
K leads to a strong decrease of the diffracted intensities 
and brings the theoretical calculations in close agreement 
with the experimental results. 

In conclusion, we computed the effective electron- 
phonon coupling constant using first-principles density 
functionnal theory and found that it strongly depends 
on the electron temperature. In addition, we reproduced 
the time evolution of Bragg peak intensities measured by 
Sciaini et aZ[llJ in FED experiments and found that the 
decay of these intensities is not only due to the increase 
of the lattice temperature but also to the redshift of all 
optical modes arising from an increase of the electron 
temperature. 

Calculations were performed using HPC resources from 
GENCI-CINES (project 095096). We aknowledge G. Sci- 
aini and R.J.D. Miller for providing us with their exper- 
imental results and also for enlightening discussions. 
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